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In this paper we show that the virial expansion up to third order for the osmotic pressure of 
a dilute polymer solution, including first-order perturbative corrections to the virial coefficients, 
allows for a scaling description of phase-separation data for polymer solutions in reduced variables. 
This scaling description provides a method to estimate the 0-temperature, where demixing occurs 
in the limit of vanishing polymer volume fraction <j> and infinite chain-length N, without explicit as- 
sumptions concerning the chain- length dependence of the critical parameters T c and (j> c . The scaling 
incorporates three limiting regimes: the Ising limit asymptotically close to the critical point of phase 
separation, the pure-solvent limit, and the tricritical limit for the polymer-rich phase asymptotically 
close to the theta point. We incorporate the effects of critical and tricritical fluctuations on the co- 
existence curve scaling by using renormalization-group methods. We present a detailed comparison 
with experimental and simulation data for coexistence-curves and compare our estimates for the 
G-temperatures of several systems with those obtained from different extrapolation schemes. 



I. INTRODUCTION 



A quantitative description of the properties of a system in the vicinity of multicritical points remains one of the 
most interesting problems in the physics of phase transitions. Closeness to a multicritical point leads to a complex 
crossover between several lines of critical points, which cannot be described in terms of a single universality class. 
As a rule, the appearance of a multicritical point on the phase transition line is related to the interaction of two or 
more order parameters . Correspondingly, the complex crossover behavior is affected by a competition between the 
diverging correlation lengths of the fluctuations of these different order parameters. The most well-known example of 
such systems with a multicritical (specifically, tricritical) point is the 3 He- 4 He mixture 

In this paper we investigate the phase separation of solutions of high-molecular-weight polymers in low-molecular- 
weight solvents in the vicinity of the 0-point, which is defined as the critical mixing point in the limit of infinite 
degree of polymerization De Gennes B, by mapping Edwards' continuous chain model onto a Euclidean field 
theory in the formal limit of zero spin components, has argued that the 0-point in the polymer-solvent system is a 
tricritical point. A tricritical point separates lines of second-order (A-line) and first-order (triple-line) transitions. The 
states above the O-temperature on the line of zero polymer volume fraction = 0, shown by the dotted line in Fig. 
[j], correspond to the critical self- avoiding- walk singularities associated with the behavior of long (N — > oo) polymer 
molecules at infinite dilution fi^j. Later on, the mapping onto a field theory was generalized to solutions of finite 
concentration and arbitrary polydispersity f^-flOl] . The scaling field hi, conjugate to the polymer order parameter 
i/> , is zero along the A-line but becomes non-zero for finite degrees of polymerization N. The second scaling field h,2 
(scalar) conjugate to ip 2 also vanishes in the limit of infinite N and zero polymer volume fraction. The correlation 
length associated with the polymer order parameter is proportional to the radius of gyration of a polymer which 
diverges in the limit of infinite degree of polymerization. Below the (tricritical) 0-point the polymer order parameter 
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exhibits a discontinuity accompanied by a phase separation and by a discontinuity in the concentration of the polymer. 
The line of critical phase-separation points shown in Fig. [I] is a nonzero-field critical ( "wing" ) line originating from 
the tricritical point. The order parameter for the fluid-fluid phase separation, associated with the volume fraction 
of monomers <p, and the polymer order parameter ip belong to different universality classes. Tricriticality emerges 
as a result of a competition between these two order parameters and exhibits mean-field behavior with logarithmic 
corrections JTT]— p3|| . Specifying the precise physical meaning of the polymer order parameter tp is a bit complicated 
(see Ref. [[| p. 287f). In analogy to the A-transition in 4 He, one can view ij> (r) as an operator for initiation or 
termination of a polymer chain at a point r, thereby relating it to the concentration of polymer endpoints. On the 
level of two-point correlations, one finds that the transverse correlations of ip are related to the correlations of the 
ends of a single polymer, while the longitudinal correlations (in the direction of hi) are related to the correlations 
between all chain ends 0^ ■ A full description of the phase separation near the tricritical point should incorporate a 
crossover between Ising critical behavior on the wing critical line and tricritical behavior close to the 8-point. 

The separation of a polymer solution into two coexisting phases has been qualitatively explained long ago by Flory 
and Huggins [fblf . However, the Flory-Huggins theory describes the phase separation only qualitatively, because it is 
based on a mean-field lattice model which neglects fluctuations. It leads to the following Helmholtz free energy F of 
mixing of polymer chains and solvent molecules |||l5| : 

/=(l-0)ln(l + In (1) 



flk B T J K T/ v N IT ' 

where f2 is the number of lattice sites in the system, k B Boltzmann's constant, T the temperature, while 4> is the 
volume fraction of the polymer, and 6 is the O-temperature. We obtain the Gibbs free energy G via the Legendre 
transformation G = F + PV, where P is the pressure and V — fll 3 , with I being the length of a lattice cell. From the 
Gibbs free energy one can calculate the chemical potential /ip of the polymer and /is of the solvent [jl5[ 

fi P (^P,T) = ^ P (cj ) ,P,T)-4(T) = k B T(f + (l-cj ) )^ I ) +P1 3 , (2) 



fi s (<f>, P, T) ee /i S (0, P,T) - M °(T) = k B T ( / - M ) + PI 3 , (3) 



which obey the relation G/Q = #p + (1 - 0)£s- In Eqs. (§) and (§) /U P (T) and ^°(T) are the standard chemical 
potentials of pure polymer and solvent, respectively. The osmotic pressure II = —fis/l 3 + P, due to the presence of 
the polymer, is the additional pressure needed to establish equilibrium with the pure solvent across a semipermeable 
membrane: 

As(O,P + n,T)=As(0,P,T). (4) 

Using Eq. (^) we find 

In the region of phase coexistence for T < 9, the values </>i and <f>2 of the volume fractions in the coexisting phases 
are found from the condition of equal chemical potentials (or osmotic pressures) in both phases: 

As(0i,P 1 T) = As (<h,P,T), (6) 
Ap(0i,P,T) = Ap(02,P,T), (7) 

where the subscript 1 denotes the solvent-rich phase and 2 the polymer-rich phase. Coexistence curves can also be 
calculated from the free energy by a common-tangent construction 



Of 



Of 



/(6O-/(0i) (8) 



which determines the two densities <fii and (f>2 of the phase-separated system. At constant pressure and temperature 
these conditions are equivalent to the equality of the chemical potentials of the polymer and of the solvent in both 
phases. According to the Flory-Huggins theory [Q4|, the dependence of the critical temperature T c and the critical 



2 



volume fraction </> c of the polymer on the degree of polymerization TV is T c = 0/(1 + 1/^/N) 2 and <fi c = 1/(1 + V~N), 
respectively. As shown by Widom flfj], for any value of the variable 

x = \VN{1 - r/T c ) , (9) 

but only if N is large and T is close to T c , the phase coexistence in the Flory-Huggins theory can be represented 
by a scaling form: the concentration difference (f>2 — <fii , where <j>\ and (f>2 are the volume fractions of polymer in the 
concentrated and dilute phases, respectively, is given by 

where the limit x — » corresponds to the critical point when T — > T c for fixed N, and the other limit x — > oo is 
approached when N — > oo, for fixed T. 



T 




line r c (/ii) 




FIG. 1. Schematic phase diagram near the tricritical O-point: (a) shown in field variables T,h\,h2; at T — Q the A line of 
self-avoiding-walk-type critical behavior ends and a surface of first-order demixing transitions with hi 7^ emerges, which is 
bound by the wing line of critical demixing points and by the triple line located at hi = 0; (b) replacing /12 by the density 
variable <f> reveals two-phase coexistence. Slices at constant chain length N, as indicated in the figure, give the familiar 
phase-coexistence curves. 



Asymptotically close to the critical point of phase-separation (the scaling variable x — > 0) the Flory-Huggins 
theory fails, as does any other mean-field theory which neglects critical fluctuations. The critical behavior of the 
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polymer solution near T c is known to belong to the three-dimensional Ising universality class. In particular, in the 
asymptotic vicinity of the phase-separation critical point the concentration difference 4>2 — 4>i as a function of the 
reduced temperature difference r = (T — T c )/T c obeys a power law: 

4> 2 - 0! ~ 2B \rf . (11) 

However, Ising-like critical behavior is revealed only when the correlation length of the critical density fluctuations is 
much larger than the radius of gyration of the polymer chain Jl7 18 1. This fact implies that the Ising region in high- 
molecular- weight polymer solutions is confined to a very narrow temperature range near the critical point. Outside 
this region, a crossover to mean-field type behavior is observed. A Flory theory renormalized by critical fluctuations 
has been developed by Povodyrev et al. pQ | . The extent of the Ising region is predicted to narrow with increasing N 
in a manner governed by the Ginzburg criterion, and disappears entirely in the limit of infinite N. For large N, the 
interplay between |r| and 1/y/N will drive the system from asymptotic Ising critical behavior, when ItIvJV 3> 1, to 
asymptotic tricritical behavior, |r| \/~N <C 1, through a region of intermediate crossover behavior. 

The approach developed in Ref. |2(| contains most essential features of real systems, namely: the crossover from 
mean-field to asymptotic Ising-like behavior, the crossover to the (tricritical) 8-behavior, and the variation of the 
effective critical exponents upon changing the degree of polymerization. However, the restrictions imposed by the free 
energy (Q) of the Flory-Huggins theory are too tight to represent actual data within experimental accuracy. This is 
why here we shall use a more general approach based on the virial expansion of the osmotic pressure n, which we shall 
renormalize so as to include the effects of fluctuations. We then try to describe existing experimental and simulation 
data on phase separation in polymer solutions with this approach. 

This paper is organized as follows. In Sec. || we briefly review the connection between polymer solutions and 
field theory and show that the Helmholtz free energy has a scaling form in the mean-field approximation and also if 



one includes contributions from the first-order perturbation theory. In Sec. Ill we describe how critical fluctuations 
renormalize the free energy and in Sec. IV we incorporate the effects of tricritical fluctuations. Sec. |v| contains 



a detailed comparison of the renormalized theory with the phase-separation data available from experiments and 



simulations. In Sec. VI we summarize our findings. 



II. MEAN-FIELD DESCRIPTION 

Utilitising des Cloizeaux's mapping between Edwards' model of polymer solutions in the case of nonzero concen- 
tration (j) and a 0(n) symmetric field theory in a magnetic field of magnitude hi in the limit n — > 0, we find the 
following relations between the polymer observables <f>, II, c p = <fi/N and the effective potential T(ip,h2) of the field 
theory gg| 

Cp= 2^' (13) 
n = -#i-r, (i4) 

where c p is the concentration of polymer molecules, N is the mean degree of polymerization. For the effective potential 
r, we start from the mean-field expression 

r = h 2 ^ 2 + \ip A + vi> 6 , (15) 

which one may interpret as a Landau expansion in terms of ip, but which also emerges as the result of a zero-loop 
("tree") approximation in a field theory or in the Edwards model |l^]. We should mention, that des Cloizeaux's 
mapping j?J demands a polydisperse grand canonical ensemble with an exponential chain-length distribution, while 
the experiments that we want to describe are performed with nearly monodisperse samples. A broad chain-length 
distribution in the experiments would seriously affect the demixing transition, since the chain-length distributions 
may differ in the phases below the demixing temperature, and thereby may even change the nature of the transition. 
Polydispersity effects have been incorporated in the theoretical framework for the excluded-volume region, where 
the effective monomer-monomer interactions are repulsive and therefore A is positive |^,^4|. Fortunately, in the 
tree approximation the effective potential T and thereby the osmotic pressure n are independent of polydispersity. 
Also the first-order correction for LI yields a change of only about 12% when changing from a monodisperse to an 
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exponential chain- length distribution p2| . Therefore, using the mean- field approximation for T, we confidently apply 
the theory to experiments on nearly monodisperse solutions. 

Minimizing Eq. (|lj) for the osmotic pressure IT with respect to ip we find 

h i = g-r= ^( 2 ^2 + 4A?A 2 + 6«V 4 ) , (16) 
for the field hi conjugate to the order parameter ip. Evaluation of (JlJ) and ([l3]) gives 

= V 2 (17) 
h 2 = jz + 2A?/> 2 + 3mp 4 . (18) 

Using ( |l7| ) and in Eq. (|l4|) we arrive at the familiar expansion of the osmotic pressure near the O-point 

n = + H 2 + 2v^ 3 • (19) 

Equation ( |l9| ) is the virial expansion of the equation of state. Note that integration of Eq. (|J) using Eq. ( |l9| ) for ft 
leads to a free energy 

/= ll n 0+A0 2 + V 3 , (20) 

similar to Flory's expression ([!]), with the (l — (f>) ln(l — <fi) term expanded and with the expansion constants replaced by 
two system-dependent interaction parameters A and v. The parameters A and v are expected to be analytic functions 
of the temperature with A changing its sign at the 0-temperature, leading to the polymer-chain-collapse transition, 
while in the tricritical scenario v is positive in order to stabilize the system at a finite density. To lowest order we can 
write 

T - 6 

A = Ao-q-, (21) 

where Ao > is a system-dependent constant, so that A is negative below the O-point and is positive above the 
O-point. The lowest-order approximation for v is simply a constant independent of temperature. 

The lambda line is reached for T > O and hi = h<z = (see Fig. |l|). At the triple line (coexistence curve in the 
limit N — > oo) for T < O we also have hy = but hi ^ 0. At the ©-point all three fields hi,hi and A vanish. In 
the semi-dilute (sd) limit of polymer concentration c p — > at fixed <p (i.e., AT — > oo) the osmotic pressure becomes 
— 4> 2 (X + 2v<f) and one easily finds from the equality of the osmotic pressure in both phases: 

<^i=0 and 0sd = 02 = -7r (22) 

zv 

for the polymer volume fractions in the coexisting pure solvent phase and semi-dilute phase. Substituting this result 
in Eq. ( |l6| ) we note that we approach the symmetry plane hi = for N — > oo. In the semi-dilute limit all polymer 
molecules aggregate in the polymer-rich phase and their overlap s — RgC p tends to infinity, because the mean radius 
of gyration R g diverges faster than c p vanishes. The critical demixing point can be found from the stability conditions 



T,N \ / T,N 



Insertion of Eq. (jnj) yields the critical volume fraction 



6v ^QvN ' 



(23) 



(24) 



^From Eqs. ( |22| ) and (|24|) we can form the universal ratio (j) s d/<t>c — 3, which, since A ~ T — O, measures the ratio of 
the slopes of the wing critical line and the phase separation boundary in the semi-dilute limit. The result </> s d/0 c = 3 
coincides with Flory theory but differs from the result s d/0c = 5/2 [ p5|]2(| that one finds from the Landau expansion 
by keeping h 2 constant instead of N in the calculation of the wing critical line. To our knowledge the only attempt 
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to independently obtain this universal ratio was made in simulations by Frauenkron and Grassberger p7| , p8[ . Their 
result s d/0c ~ 2.9 — 3.2 supports the validity of a calculation at constant N. 
Using the scaled variables 

A0 = tlh- , (25) 

0c 

AA=^, (26) 
and the mean-field result (|24|), we find that the osmotic pressure ( |l9| ) obeys the scaling form 

An = n 7 Hc = -3AA(1 + A0) 2 + (A0) 3 , (27) 

with il c = 2v0 3 and without explicit system-dependent parameters. To obtain the scaling equation for the free energy 
we subtract its regular part 



df 

/reg = f{(f>c) + TT7 



(0 " 0c) , (28) 



which neither does contribute to the susceptibility nor affects the calculation of coexistence curves. With this sub- 
traction we find 

J~f re &—orr-i , a!m„m i a 2\ a 1 /a\, 1 WaI\2 , 1 r a l\3i 



A/ = = 3[(1 + A0) ln(l + A0) - A0 - I AA + - I (A0) 2 + -(A0) 3 ], (29) 



and the coexistence curve is calculated from the conditions 



dAf 



dAf 



,i <9A0 



dA<j) 

The Ansatz (|l|) for the temperature dependence of A gives 



A/(A0 2 ) - A/(A0!) _ 



A ^ 2 A0 2 -A0! 



T — T 

AA = ^, (31) 

for the scaling variable AA. All coexistence curves are expected to collapse onto a single scaling curve, when plotted in 
terms of the variables (T — T C )/(T C — Q) and (0— C )/0 C . Such scaling behavior was already proposed earlier by Izumi 
and Miyake |37[ based on homogeneity arguments, but the quality of the scaling was not very good. As mentioned 
earlier j25), the scaling variable ( |3l| ) for critical temperatures T c close to the 8-point is rather sensitive to the value 
of the 0-temperature. The extrapolation of = limjv-^co T C (N) is notoriously difficult due to poorly controlled finite 
chain-length effects and effects of polydispersity. This is why in the earlier work |25| both — C and T — T c were 
scaled by C . In practice, it means that T c — was replaced by an empirical function of C taken from experiment. 
Moreover a second-order term (quadratic in T — T c ) was added to account for nonasymptotic effects at lower degrees 
of polymerization. Such an approach yields an almost perfect description of the phase separation but contains too 
many empirical features. In this work we try to avoid any empirical assumptions and keep T c — as a scaling factor 
while adjusting the value of the 0-temperature. To derive the scaling we use the zero-loop approximation a 2 = A and 
03 = 2v for the second and third virial coefficient in the virial expansion ([l^) of the osmotic pressure. Perturbative 
corrections in the two- and three-point couplings A and v have been calculated up to two- loop order [Q, but since 
this leads to divergent series for the virial coefficients, which one does not know how to resum, we restrict ourselves 
to incorporating only the first-order perturbation as calculated by Duplantier The calculation to first order in 
the coupling constants A and v modifies only the second virial coefficient, which for a monodisperse solution reads in 
our notation 

24i> 

W2 = A ~ (27r)3/ 2 iW2 ' (32) 

Solving the critical point conditions (^3|), we find that the scaled Eq. (^) for the osmotic pressure still holds, but 
with a modified scaling variable 



G 



AA = f i _ 4 ^ = Lzll L 4 , f33 ^ 

A c V 1 (2^)3/20 c 7Vi/2; T C -Q\ L (271)3/2^1/^' 



Inclusion of the effects of tricritical fluctuations, to be described in Sec. [V, introduces a nonuniversal parameter 
p{v^) in the denominator of the correction term in Eq (|33|). One can also argue, that such a factor should already be 
accounted for in the unrenormalized result, since the microstructure of the real polymer differs from that of the model 
underlying the calculation. Thus the relation between the chain-length of the actual polymer and the chain-length 
used in the calculations contains a nonuniversal factor. This consideration leaves us with two fitting parameters 
and c = 4/(27Tyo(?;pJ) 3 / 2 , and with the chain-length N — M w /M m measured as the number of monomers in a chain 
(M m being the monomer molecular weight). 

III. CRITICAL RENORMALIZATION 

To incorporate the influence of long-range critical fluctuations, which strongly affect the shape of the coexistence 
curve in the vicinity of the critical demixing point, we employ a procedure developed by Chen et al. pl|,p2"| using 
renormalization-group matching [^9| to implement the crossover between mean field and critical behavior. The free 
energy (^9|) can be expanded around the critical point, yielding 

A/ = -3AA(A0) 2 + i(A0) 4 + .... (34) 

The renormalization-group theory for the scalar tp 4 field theor y |30| provides a proper tool to account for critical 
fluctuations that give rise to scale invariance at the critical point |31[| . The scaled free-energy of a system with order 
parameter M and reduced temperature t in mean-field approximation can be expanded as 

A/=^M 2 + ^M 4 + ..., (35) 



where the four-point coupling u is normalized by its fixed point value u* = 0.472 [g0|,^2| and where A is a cutoff 

length-scale. Since all system-dependent parameters are scaled away in Eq. (|35]), we allow for two amplitudes c p and 
c t in the relations 

t = c t AX, M = CpA(j), (36) 

connecting our physical variables to those of the field theory. Equating the expansion coefficients in Eqs. (|34|) and 
(|35|) we find the relations 



etc 2 = 6 , u*uA = — , (37) 
c p 

restricting the four parameters u, A, c t , c p to only two, which can be varied independently. The above expansions do 
not depend on the degree of polymerization. It is interesting to compare expansions ( |34"| ) and (|35|) with conventional 
Landau expansions of the unsealed free energy p3] 

/ = ia o ^^(A0) 2 + i Uo (A0) 4 + ... (38) 

and with that of the reduced variables 

1 , %9 u*uA s , , 

f=-t{Mf + —u {Mf + ..., (39) 

where 

t = c t 1 -^-, Af = c p A0, (40) 

J c 

and 

c t c 2 p = a , u*uA = ^ . (41) 
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In the Flory model we find in the large N limit ao = 2T c /0Ao — > l,iio = 12w/0 c — * 2/0 c , hence Ao = 2 and 
V = 1/6. Furthermore, since c t = c t o/VN,A = Aq/*/N, and c p = Cp^/N 1 / 4 JlS| ], comparison with Eqs. (|^) and ( |37| ) 
yields: A = A, c^o = c p , and c f0 = cj/6. Note, that these parameters, as well as u, do not depend on the degree of 
polymerization. 

Renormalization proceeds by replacing the physical variables in Eq. (p5[) by renormalized ones defined as 



t^t x =tT, M -> M x = MX> 1/2 , u -> u x = , (42) 

where the rescaling functions T, 2? are integrated renormalization-group (RG) exponent functions and hi is the 
integrated Wilson flow function pl[ . We can approximate these functions with good accuracy in terms of a crossover 
variable Y by §§] 

r = y (2-i/i/)/« j T> = Y -*/ u , U = Y 1/uJ : (43) 

with cj = A B /i/, where 77 = 2 - 7/1/ = 0.033 ± 0.003, 7 = 1.239 ± 0.002, zv = 0.630 ± 0.001 , and A s = 0.51 ± 0.02 
^o| , ^3| -[36l are the universal critical exponents of the asymptotic power laws for the correlation function, susceptibility, 
correlation length, and for the Wegner correction to asymptotic scaling, respectively. For the free energy itself, it is 

necessary to perform an additional additive renormalization by adding the "kernel" term — — P/C with 

K = -^(y-"A^ - l) (44) 
auk V / 

to the free energy density. In Eq. (f44| ) a = 2 — = 0.110 ± 0.001 |3^,^ |35| is the universal critical exponent for the 
heat capacity. The crossover function Y is defined implicitly by the equation 

1 - (1 - u)Y = u[l + (A/k) 2 ] 1 ^ 1 ^. (45) 

which evaluates the integrated flow equation for the running coupling constant u(l) at a specific matching point I = I* 
where one recovers the mean- field expression for the free energy p2) . The parameter k is defined as 

~ K >= 92 f- , (46) 

where A/ x = A/(t x , M x ). It measures the distance from the critical point. Note that the unsealed k = k/N 1 / 4 depends 
on the degree of polymerization, being inversely proportional to the correlation length £ Jl9|,0 ■ Asymptotically close 
to T c ,k 2 — * and Y ~ (k/wA) u , yielding the Ising asymptotic behavior fl^j . Close to the 0-point k 2 ~ AA = 
(T — T C )/(T C — O), diverging as T c — > Q and driving the crossover function Y to its mean-field limit Y — 1. In 
unsealed variables, the crossover to mean- field 0-point tricriticality is driven by A ~ iV -1 / 2 , assumed to be inversely 
proportional to the radius of gyration. To obtain the crossover to mean-field tricriticality, instead of expanding the 
free energy A/, we now renormalize the full expression ( p9| ) to keep the proper low-density limit away from the 
Ising critical point, which is contained in this expression. A similar approach was used earlier to incorporate critical 
fluctuations into the Van der Waals equation |59| and into the Flory-Huggins model ptJ . Since Eq. (|2^) contains only 
AA and Acj) as variables, corresponding to t and M, but no explicit coupling u, we incorporate the renormalization 
factor U into the factors renormalizing AA and A</> in a way, that reproduces the renormalization of ([55]) when j2S| ) 
is expanded. This leads to the renormalization transformations 

AA -► AA X = TU~ 1/2 A\ = SAX , A0 A0 X = V 1/2 U 1/4 A4> = TZA4> , (47) 

and an accordingly modified definition of the crossover parameter k 

k 2 =U^x x \ (48) 
c 

l> 

where \~ 4 is the renormalized ("crossover") inverse susceptibility 

„_! 9 2 A/ , aX Al , J 1 nfAX , 1 



X* = ^77lW( AA - = 3 \ -Zi - 2 ( AA * + 9) + A ^ ■ ( 49 ) 

oyAcpy \1 + A0 X * 
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We observe that Eq. ( fig ) can be solved for A</> x and after eliminating k 2 with Eq. fllq), we obtain A0 X as a function 
of the crossover variable Y and AA. We then numerically solve the two-phase coexistence conditions 



dA<p 



a i dA(f> 



A/ x (A0 x2 )-A/ x (A0 xl ) 
A02 - A0i 



for Yi and Y2 at a given temperature and thereby find the reduced densities A0i and A02 of the coexisting phases. 
The procedure outlined above implements the crossover between the critical regime and the tricritical regime. In the 
critical regime the correlation length of concentration fluctuations is larger than the radius of gyration of a polymer, 
which at low enough concentration sets the scale for tricritical fluctuations. Thus near the critical demixing point 
the coexistence-curve has an Ising shape. In the tricritical regime the radius of gyration is larger than the correlation 
length of concentrationl fluctuations and the coexistence curve becomes triangle shaped. One can quantify the location 
of the crossover region by a scaled Ginzburg number, which is given in terms of our crossover parameters by 

iVc= 5 o^, (51) 

with go = 0.031 |po| . Fitting the theory to the experimental data as described in Sec. [v], we find (in agreement 
with light-scattering experiments 0) u = 1 and c p = 1.6, and with the use of Eq. ( |37| ) we obtain Nq = 0.05. For 
AA <C Nq we find Ising-type behavior and for AA ^ Nq we find tricritical behavior. In unsealed variables, the actual 
Ginzburg number is Nq — Nq\ c , which in the Flory model becomes Nq = Nq/^/N J|(|. 

IV. TRICRITICAL RENORMALIZATION 

Renormalization-group treatment of the n-vector model in the limit n — > 0, describing tricritical behavior , reveals 
that the upper critical dimension for the tricritical point is d c = 3 f40]] . Thus, in dimension d = 3 we expect mean-field 
behavior at the tricritical point, implying that O-point polymers on a large scale behave effectively like random walks 
with a mean-field correlation-length exponent v = 1/2. In the vicinity of the tricritical point one finds logarithmic 



corrections to mean-field behavior. While for the scaled free energy (29) tricritical renormalization factors cancel to 
a large extent in the scaled variables, the chain-length dependence of the critical parameters </> c (iV) and T C (N) is 
modified by tricritical fluctuations. To lowest order the tricritical renormalization-group mapping reads |l2| 

A -> A R = Aa«)4 /n , (52) 

v ~+ VR = srbv - (53) 

0-0a = M«a)(l + 0(vR)), (54) 
N^N K = Np(vi)(l + 0(v R )), (55) 

where <j(v^) and p(fp) are nonuniversal integration constants depending on the start value of the renormalized 
three-body interaction. Renormalizing Eq. ( p4| ) one finds the asymptotic relations 

N © \fN 



As was shown by Hager and Schafer 113], the region where the asymptotic relation (|53| ) for wr, is valid is limited to 
chain lengths way beyond of what is within the reach of present experiments or simulations. Despite recent progress 
ijllf , the tricritical Wilson flow function in the crossover region is still known with much lower accuracy as compared to 
the critical case. We thus choose to eliminate the running coupling v^{N) in the critical-point conditions, so as to be 
able to test the tricritical predictions without using some approximation for wr. Without tricritical renormalization 
we find from Eqs. ( p^ ) and (|Tj) the mean-field relation 

^-^ ~ — (57) 

e ^n ■ (57) 

If we first renormalize Eq. (|2~i| ) using the mapping (j|^)-((5|), we find 
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e " ~ ^nt/h ' (58) 

as the renormalized mean- field result, and if we include the first-order correction of Eq. (^) with A — > A in accordance 
with Eq. p2|, we obtain the relation 



e-r c 1 



1 - 7^-7-oTTT77-T7T77 ■ (59) 



Besides the O-temperature and the parameter p(v#), these relations contain only directly measurable quantities and 
we shall test their validity in the next section. Let us now obtain the chain-length dependent Ginzburg number 
Ng(N) — A/gA c - With the use of the asymptotic tricritical result ( |56| ) we find 

In" 3 / 11 N 

N G (N) = Nq\ c j=— . (60) 

v TV 



Thus the Ising region shrinks to zero, when we approach the 0-point (N — > oo) but the classical A -1 / 2 dependence is 
now modified by tricritical fluctuations. By comparing the expansion coefficients of the unsealed free energy / — / rcg 
to Eq. (|3^) and with the use of Eq. ( |5l| ) and the asymptotic tricritical results ( |56] ) , we find the chain-length dependent 
parameters 



1 /3 

fi(W )A( W > = ^pf&) (61) 

A c y .go J VN 

c p (N) = fcW 8 (^-) V6 ~ v^ln- 1 / 12 A . (63) 

Apart from the logarithmic terms, due to the tricritical renormalization, the chain-length dependence is the same as 
that for the Flory model po|. 



V. COMPARISON WITH EXPERIMENAL AND SIMULATION DATA 



A. Experimental coexistence curves 

We first analyze three sets of experimental coexistence-curve data available in the literature |42]-|47j| . The original 
data are displayed in Figs. 2-4. The first set of data is on the system polystyrene (PS) in cyclohexane [^-Q, with 
coexistence curves for three different molecular weights M w = 110,000,200,000,1,560,000. The next set contains 
data for seven different chain lengths for the system PS in methylcyclohexane |4q| with molecular weights ranging 
from 10,200— 719,000. Another data set covers six chain lengths for the system polymethylmethacrylate (PMMA) 
in 3-octanone |If|[l7), with molecular weights ranging from 26, 900 — 596, 000. In all cases the polymer fractions were 
reasonably monodisperse with polydispersity indices in the range M w /M n = 1.03 — 1.11, where M w and M n are the 
weight averaged and number averaged molecular weight, respectively. For the PS data we used all data points as 
published. For the PMMA data we used the critical point volume fraction <fi c = {(j>i + <j>2)/2 with cf>i and <f>2 being the 
pair of coexisting densities closest to the critical point. This leads to critical volume fractions being on the average 
4% smaller than the published ones. As a first step we fitted the parameters and p(v^), contained in the scaling 
variable AA given by Eq. (|33|), to find an optimal collapse for the data within each data set when, plotted in terms 
of the scaling variables AA and A(f>. It turns out, that p(Vjf) = 1 is a reasonable choice for all three data sets. This 
leads to the coefficient 4/(27r) 3 / 2 = 4/[27rp(vpJ] 3 / 2 = 0.254 in the correction term of Eq. ( |3^ ) and, consequently, in 
Eq. d59|). The contribution from the correction term in Eq. (|3^) is about 13% or smaller for all data sets. The fitted 
0-temperatures are: = 309 K for PS in cyclohexane, = 345 K for PS in methylcyclohexane and = 345.65 
K for PMMA in 3-octanone. Our value for PMMA in 3-octanone lies between the values = 345.15 K |47| and 
= 346.85 K |Q, obtained by extrapolating 1/T C (N), taking A -1 / 2 corrections into account. Our estimates for PS 
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in cyclohexane and in methylcyclohexane lie slightly above the values 9 = 307.15 K and 8 = 341.95 K obtained by 
Izumi and Miyake |37|, who used the mean-field scaling variable ([bT|). The difference is due to the correction term in 
Eq. which leads to a better data collapse than in Ref. p7fl, but for slightly higher 8-temperatures. 




FIG. 2. Coexistence-curve data for PS in cyclohexane [Jt2j-|44| are shown together with the result of the renormalized crossover 
theory (full curves). The dotted curve represents Eq. with O = 307.25 K. The critical demixing points and the coexis- 

tence-curve scaling value of = 309 K are denoted by stars. An estimate for the limiting phase-separation boundary with 
4>ad/4>c = 3 is shown by the dashed line. 




FIG. 3. Coexistence-curve data for PS in methylcyclohexane ]45j are shown together with the result of the renormalized 
crossover theory (full curves). The dotted curve represents Eq. (p9|) with O = 342.75 K. The critical demixing points and the 
coexistence-curve scaling value of = 345 K are denoted by stars. An estimate for the limiting phase-separation boundary 
with 4> s &/4>c = 3 is shown by the dashed line. 
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FIG. 4. Coexistence-curve data for PMMA in 3-octanone |4^,^7|] are shown together with the result of the renormalized 
crossover theory. The dotted curve represents Eq. (|H^) with 6 = 343.15 K. The critical demixing points and the coexis- 
tence-curve scaling value of = 345.65 K are denoted by stars. An estimate for the limiting phase-separation boundary with 
4>ad/4>c = 3 is shown by the dashed line. 



Compared to other methods, that use mean-field extrapolations of the critical parameters 4> c and T c leading to 
values = 305.6 - 308.4 K (see Ref. 0| p. 163) for the system PS in cyclohexane, our value is at the upper end 
of the estimated range. This method of fitting the O-temperature has the advantage, that no explicit TV-dependence 
for T C (N) needs to be assumed. The scaled data are displayed on a linear scale in Fig. || and in a logarithmic scale 
in Fig. U Note that the one data set of Ref. [[l6| with M w = 26, 900, which is significantly off the scaling curve (see 
log plot), has also the largest polydispersity index M w /M n = 1.11 of all samples. It is interesting to observe that the 
data for all three systems with good accuracy collapse onto a single scaling curve, as is suggested by the solution of 
the unrenormalized coexistence conditions ([50]) indicated by a dashed line in Figs. || and [|. This is not necessarily 
to be expected, since the critical renormalization introduces two additional parameters u and c p which may have 
different values for different solutes and solvents. One clearly sees, that the unrenormalized scaled coexistence curve 
fails to reproduce the proper Ising-type singularity (|ll]) with (3 — 0.325 close to the critical point, but instead has 
the mean-field exponent (3 — 1/2. In our second step in fitting the data we remedy this deficiency by applying the 



renormalization procedure as outlined in Sec. III. Since the parameter u was already found to be close to unity in a 
recent evaluation of light-scattering data above T c [l8| , and its variation does not affect the coexistence curves very 
much, we fixed it to u = 1, leaving c p as the only fit parameter of the crossover theory. The result of our fit with 
c p = 1.6 is displayed in Figs. 2-4 in terms of unsealed variables and in Figs. 5 and 6 in terms of scaled variables. Our 
scaled crossover formulation nicely reproduces the Ising singularity, but shows some deviations in the crossover region 
which cannot be removed by tuning c p . One can think of a plethora of possible sources for such a deviation, since in 
our approach we neglected higher-order terms in expansions at several stages in our calculations. We truncated the 
virial expansion at order 3 , included only first-order perturbative corrections in A and v, and used the lowest-order 
approximations for the temperature dependence of A and v. The influence of higher-order terms in the expansion of 
the temperature dependence of A and v and similarly of higher-order terms in the perturbation theory causes only a 
change of the scaling variable AA, that cannot account for the deviation. Note that both coexisting densities 4>i and 
c/>2 of the renormalized fit are shifted to higher densities as compared to the experimental data. Higher-order terms 
in the virial expansion, like a negative <fi -term can induce a shift of both coexisting densities in the double-tangent 
construction (50) to lower values. Since we use the theory to volume fractions up to < 0.4, it is likely that such 
higher-order terms in the virial expansion contribute to the observed deviations. 



12 




( * - i ) / <> c 

FIG. 5. Linear scaling plot of phase-coexistence curves for PS in cyclohexane p^-[l4|] (full symbols), PS in methylcyclohexane 
(open symbols) and PMMA in 3-octanone |46t[l7| (open crossed symbols). The dashed curve represents the values calculated 
from the unrenormalized scaled free energy (129) and the full line represents the result of renormalized crossover theory with 
u — 1 and c p = 1.6. 
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(T - T )/(T - 0) (1 - clifl N 1/2 ) 

FIG. 6. Logarithmic scaling plot of phase-coexistence curves for PS in cyclohexane p2|-|44[ (full symbols), PS in methylcyclo- 
hexane (open symbols) and PMMA in 3-octanone |l 6| |47 | (open crossed symbols). The dashed curve represents the values 
calculated from the unrenormalized scaled free energy (G3) and the full line represents the result of renormalized crossover 
theory with u = 1 and c p — 1.6. 
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B. Scaling of the critical parameters 



In the previous section we assigned values for the 0-temperature that optimized scaling of the coexistence curves. 
An alternative procedure to estimate the O-temperatures is by extrapolating N — ► oo in accordance with Eqs. ( p38| ) 
and (||). The values thus obtained for (6 - T c )/Q are plotted in Fig. | as a function of C T 3/11 TV" 7 / 11 . The solid 
lines indicate the asymptotic linear relations in accordance with Eq. (|58|). The slope of this linear relation was found 
to be 2 for PS in cyclohexane, 2.5 for PS in methylcyclohexane and 3.5 for PMMA in 3-octanone. The resulting 
estimates for the O-temperatures are O = 307.25 K for PS in cyclohexane, © = 342.75 K for PS in methylcyclohexane 
and O = 343.15 K for PMMA in 3-octanone. They all lie about 0.7% below the values found from coexistence-curve 
scaling. Those for polystyrene compare well with other extrapolations of critical-point data p^j . We believe that the 
origin of the discrepancy again lies in the neglect of higher-order terms in our calculations. To discriminate between 
the relations ( |57j ) - ( p9| ) we form a ratio A, by dividing the right-hand side of each equation by its left-hand side. The 
resulting ratio should be constant. The values of the ratio A for Eqs. (|57|) - (|5^) and the three data sets are displayed 
in Fig. ^, with the left-most values being normalized to unity within each set and with an additive offset to separate 
the data sets. One clearly sees, that the unrenormalized prediction of Eq. ( |57j ) gives the worst fit and this cannot 
be improved by lowering the 0-temperature, since then the curves start bending upwards close to 0-temperature. 
The renormalized fits for PS in cyclohexane and PMMA in 3-octanone are indeed constant within the scattering of 
the data, but the most precise data on PS in methylcyclohexane clearly show a residual slope which can be easily 
accounted for by introducing a correction (0 — T c ) 2 . This, together with the mismatch to the ©-temperatures from 
coexistence-curve scaling, indicates the presence of second-order terms. Our fits for the critical-point scaling are 
included as dotted lines in Figs. &0. We also included the O-temperatures obtained from coexistence-curve scaling 
and a guess for the limiting phase boundary </> s( j, assuming the ratio <f> s d/ 4>c — 3. Note that a value of 4'sd/4>c — 5/2 is 
not in accord with the experimental data unless one assumes the presence of large corrections to the asymptotic value. 
We think this is unlikely since the value of the ratio is not altered by the leading order of tricritical renormalization 





0,07 



FIG. 7. The reduced temperature difference between the critical point and the 0-point is plotted versus C 3//11 7V~ 7/ ^ 11 . 
The symbols denote the critical parameter data for three experimental systems. The lines indicate the asymptotic linear 
approximation in accordance with Eq. (|58[) . 
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FIG. 8. The ratio A (right side / left side of Eqs. plotted against the reduced difference between the critical 

temperature and the 0-temperature. The experimental data for three different systems are separated by an offset of magnitude 
1. 
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FIG. 9. Linear scaling plot of simulated phase-coexistence curves ^] for self-avoiding walks with nearest-neighbor interaction 
on a simple cubic lattice. The dotted curve represents the values calculated from the unrenormalized scaled free energy ( |2s| ) 
and the solid curve represents the renormalized crossover theory with u = 1 and c p = 1.5. 



C. Simulations 



We have also analyzed a set of coexistence-curve data obtained from computer simulations of self-avoiding walks 
with attractive nearest-neighbor interaction on a simple cubic lattice |49[ along the same lines as was done for the 
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experimental coexistence-curve data in Sec. VA. We considered seven sets of coexistence-curve data with chain 
lengths ranging from 64 to 1, 000 monomers. Figure || shows a linear scaling plot of the simulation data and Fig. [h] 
a logarithmic one. The optimal data collapse was obtained for O = 3.754, measured in units of the nearest-neighbor 
interaction energy, and c = 0.51. The value for O again is 1% higher than the best value O = 3.717(2) obtained from 
extrapolation of 1/T C [^Tj. In general the quality of the data collapse is less satisfactory, compared to the experimental 
data. Especially on the semi-dilute branch the data for N = 800 and N = 1, 000 deviate significantly from those for 
smaller chain lengths. One main difference between experiment and simulations is that the simulations were done 
in rather small boxes. These boxes were chosen to be at least four times larger than the maximum value of the 
radius of gyration. However, with the mean squared end-end distance of a random walk being a factor -\/6 larger, 
this still implies that a single polymer stretches over a significant portion of the box, causing finite-size effects. Also 
the rather small and varying total number of polymers, ranging between 80 and 350 for different simulations may 
cause problems in obtaining reliable bulk estimates, since large surface contributions are to be expected. On the 
other hand, we observe a clear breakdown of coexistence-curve scaling for the simulation data for the shorter chain 
lengths 8, 16 and 32. This is to be expected, since we then reach volume fractions where our virial expansion is no 
longer accurate. Furthermore, the symmetric coexistence curve of the Ising model for N — 1 cannot be scaled onto 
a asymmetric scaling curve with a finite reparametrization. The two fit parameters u and c p of our crossover theory 
are not sufficient to provide a decent fit to the scaling curve, probably due to the reasons already discussed for the fit 
of the experimental data. More work has to be invested in order to decide, whether the observed mismatch is due to 
a limitation of the coexistence-curve scaling or can be attributed to systematic deviations in the simulation data. 




1 1 1 1 — 1 — 1 1 1 1 i 1 1 1 — 1 — 1 1 1 1 i 

0,01 0,1 1 

(T - T)/(T - 0) (1 - difl N 1/2 ) 

FIG. 10. Logarithmic scaling plot of simulated phase-coexistence curves [^9) for self-avoiding walks with nearest-neighbor 
interaction on a simple cubic lattice. The dotted curve represents the values calculated from the unrenormalized scaled free 
energy (EST) and the solid curve represents the renormalized crossover theory with u — 1 and c p = 1.5. 



VI. CONCLUSION 



We find that the scaling description for coexistence curves of polymer solutions, including first-order perturbative 
corrections, provides a decent collapse of the available experimental data onto a single scaling curve at fitted 0- 
temperatures which are systematically about 0.7% higher than those obtained from an extrapolation of critical point 
data to infinite chain-length. Incorporation of the critical fluctuations via crossover RG theory leads to a decent 
description of the scaled coexistence curve by fitting essentially only one parameter c p = 1.6. We attribute the 
remaining deviations in the fit and the discrepancies in the 0-temperature estimates mainly to two approximations in 
our approach. First, we truncated the virial expansion after the (j> 3 term, neglecting higher-order terms that contribute 
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at larger volume fractions, and will perturb the scaling description as they do in the Flory-Huggins model. Another 
source of corrections are higher-order terms in the temperature dependence of the bare parameters and higher-order 
corrections of perturbation theory, which have similar effects on the scaling variable AA. We find that both simulations 
and experiments are pointing towards an universal ratio 4> s d/4'c = 3, implying that <f) c has to be calculated at constant 
chain length N rather than at constant field /12, which would lead to 4> s d/4'e — 5/2 in the Landau theory. The scaling 
description is less satisfactory for coexistence-curve data provided by simulations. Concerning the 0-temperature 
extrapolations wc find from the simulation data the same pattern as for the experimental data, with the value of O, 
obtained from the scaling description, being 1% higher than the best value obtained from a extrapolation of l/T c (N). 
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